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We study domain coarsening of two dimensional stripe patterns by numerically solving the Swift- 
Hohenberg model of Rayleigh-Benard convection. Near the bifurcation threshold, the evolution of 
disordered configurations is dominated by grain boundary motion through a background of largely 
immobile curved stripes. A numerical study of the distribution of local stripe curvatures, of the 
structure factor of the order parameter, and a finite size scaling analysis of the grain boundary 
perimeter, suggest that the linear scale of the structure grows as a power law of time t 1 ^ , with 
z = 3. We interpret theoretically the exponent 2 = 3 from the law of grain boundary motion. 
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Equilibrium layered phases (characterized by a uni- 
form wavevector ko ^ 0) are often found in systems 
with competing short and long ranged interactions |Q. 
Related structures, commonly referred to as stripe pat- 
terns, also appear in systems driven outside of thermody- 
namic equilibrium (e.g., Rayleigh-Benard convection or 
parametric surface waves near onset &), After chang- 
ing rapidly a control parameter across a transition or 
bifurcation point, a uniform state become unstable and 
configurations with locally ordered stripes appear. Given 
the underlying translational and rotational invariances of 
the system, spontaneous evolution leads to a macroscopic 
sample comprising a large number of grains or domains, 
each relatively uniform, but oriented along an arbitrary 
direction, as well as to a large density of defects such as 
grain boundaries, disclinations and dislocations. Under- 
standing how this structure orders with time, and how 
the motion of interacting defects contributes to the coars- 
ening rate is the main focus of this paper. 

Numerical studies of model equations in two spatial 
dimensions j3|-0, as well as recent experiments involv- 
ing thin films of block copolymers |sj support the idea 
that the time evolution of layered phases after a quench 
is statistically self-similar (the statistical self-similarity 
hypothesis asserts that after a possible transient, con- 
secutive configurations of the coarsening structure are 
geometrically similar in a statistical sense). As a conse- 
quence, any linear scale of the structure (e.g., the average 
size of a domain or grain) is expected to grow as a power 
law of time R(t) ~ i 1/z . 

Coarsening of layered phases is not yet well under- 
stood. On symmetry grounds, layered phases can be 
classified as smectics ||[|. Hence, by analogy with coars- 
ening studies of nematics (hJ and 0(N) vector models 
with a nonconserved order parameter [nTl , one would ar- 
gue that self-similar coarsening is to be expected with 
z = 2. Although, the possibility of a long time cross over 
to z = 2 has in fact been considered [0, numerical evi- 
dence has consistently pointed at values of z in the range 



z = 4 — 5 [^-[7|. More importantly, the self-similarity 
hypothesis itself has been questioned as different linear 
scales yield different values of z HQ. Furthermore, and 
in contrast with related research on well understood sys- 
tems that order at ko — [p^|Jl^ 1, the value of z appears 
to be modified by the presence of thermal noise. 

We present here a numerical investigation of domain 
coarsening for the Swift-Hohenberg model of Rayleigh- 
Benard convection H|. Briefly, we focus on the region 
close to onset (e — * 0, where e is the reduced control 
parameter)) and find that coarsening proceeds in a self- 
similar manner. We analyze several different character- 
istic length scales, and find that they are asymptotically 
proportional to each other. We also find that z ~ 3, in- 
dependent of thermal noise. We interpret this value of 
z from the law of grain boundary motion given in [ [L5| |. 
Further from onset (e ~ 0.25), we recover the results of 
previous studies, which we interpret as arising from non- 
adiabatic effects that lead to defect pinning. This fact 
accounts for both the slower growth seen previously, and 
its dependence on fluctuations. 

The Swift-Hohenberg model of convection in dimen- 
sionless units is \A, 
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where is an order parameter related to the vertical fluid 
velocity at the mid plane of a Rayleigh-Benard convec- 
tion cell, e is the reduced Rayleigh number, fco is the roll 
wavenumber, and £o is a constant that depends on the 
boundary conditions at the top and bottom plates. For 
our purposes, we note that £o oc 1/ko- The same model 
has been used to analyze coarsening of lamellar phases 
in a diblock copolymer M, and is otherwise believed to 
be a generic model of the kinetics of stripe formation. 
There, ip is the concentration difference between the two 
monomers jl6],[l7| . For < £ « 1, the stationary solution 
of Eq. (Q) is well approximated by a sinusoidal function 
of wavenumber ko. The transient evolution and domain 
coarsening is investigated by numerically integrating Eq. 



f 



([j]) from random initial conditions. All calculations are 
performed very close to onset (e = 0.04). Details of the 
numerical algorithm can be found in ref. jl5| . 




FIG. 1. Order parameter tj) shown in grey scale at time 
t = 15000. Eq. (|l|) is discretized on a square grid of mesh size 
Ax = 1 with 512 nodes. The wavelength Ao = 2w/ko = 8Ax. 
The reduced Rayleigh number is e = 0.04. The initial condi- 
tion has (ip) = and (ip 2 ) = 0.04. 

Figure [l] shows a typical transient configuration. The 
configuration contains a large amount of grain boundaries 
that separate domains of different orientation, as well as 
defects (such as +1/2 disclinations and dislocations). 

We first present our numerical results for several mea- 
sures of the linear scale of the structure, including the 
distribution of stripe curvatures, the order parameter 
structure factor, and the grain boundary perimeter. Fol- 
lowing initial transients, the stripe curvature k = |V • n| 
is a slowly varying field, where n denotes the unitary 
vector normal to the lines of constant ip. We com- 
pute the probability distribution function of stripe cur- 
vatures P(k, t) by considering only the subset of points 
where stripe orientation can be properly defined, i.e. the 
points that are not in the immediate vicinity of any grain 
boundary nor other defects. Far enough from a defect, 
ij)(r) = A(f) cos(k(f) ■ r + <f>), with A a slowly varying 
amplitude. By defining ((f) — if 2 + (Vip) 2 /k 2 ., one has 
((f) ~ A 2 . Note that for stationary parallel stripes of 
wavenumber k , ((f) — ( = 4e/3 [||. We now define 
defect free regions as those that satisfy r m < (/(o < r M, 
with r m = 0.95, rja = 1-05 (filter a) or r m — 0.97, 
tm = 1.10 (filter b). We have numerically verified that 
the values of ( corresponding to a set of moderately 
curved stripes along their transverse direction remain 
completely within the intervals defined by both filters 
a and b. By contrast, most values of ( in the vicinity of 



a grain boundary are lower than 0.90. 
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FIG. 2. Probability distribution function of curvatures 
P(n,t) averaged over 35 independent initial conditions for 
times ranging from t = 960 to 1.6 10 4 . 
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FIG. 3. Moments m n (t) of P(n,t) obtained with filter a. 
The values of n shown are 1/2,1,2,3, from top to bottom. The 
straight lines have a slope of 0.32. 

Figure || shows P(k, t) at different times, after averag- 
ing over 35 independent initial conditions on a square grid 
of size 1024 x 1024 (16 grid nodes per wavelength Ao), at 
e = 0.04 and using filter a. To check for self-similarity we 
study whether P(k, t) = t l • lz : p(nt^ • /z ) . In order to do this, 
we compute its moments m n (t) — f* ' dn K n P(n,t) 

with n c (t) defined as f$ dn P(n,t) = jK c (t)P(0,t), 
and 7 an arbitrary constant, < 7 < 1. We find that 
(77i„/r7io(t)) -1 /™ oc t 1 /* , with the value of z independent 
of n and 7, thus lending support to the self-similarity 
hypothesis. Figure § shows the results for a few values 
of n and two values of 7. In addition, the best fit to the 
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curves yields 1/z ~ 0.32 with filter a, and 1/z 
with filter b (not shown). 
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FIG. 4. Finite size scaling analysis of the total grain bound- 
ary length, with systems of aspect ratio r\ = 32, 42.66 and 64. 

We next present a finite size scaling analysis to inde- 
pendently determine the value of z |lq| . Let N4 be the 
number of grid points for which Q > Tm or £ < r m . 
The probability of a point belonging to a defect is pd — 
NdAx 2 1 L 2 , with L the system linear extent. We define a 
dimensionless defect (i.e. grain boundary) perimeter as 
C* = T] 2 pd, where rj — L/Xq is the aspect ratio. For short 
times finite size effects are expected to be negligible, and 
Pd ~ t~ x / z . We now introduce a finite size scaling ansatz, 
valid for any time t, 



£*(r ! ,t) = r ! g(t/r 1 z ), 



(2) 



with g(x) ~ x~ x l z for x <C 1. At fixed e = 0.04, we have 
numerically computed C* (t) with the help of filter a for 
different aspect ratios: 77 = 32, 42.66 and 64 (i.e. r\ — 
258/8,256/6 and 512/8; Ax = 1 in all cases). We have 
averaged the results over 500, 300 and 100 independent 
initial conditions respectively. Figure [| shows our results 
for the universal curve g(x). The value z = 3 leads to 
curves C* jr\ as a function of t/rf do not depend on rj. 

We have also analyzed the Fourier transform of the 
two-point correlation function of the order parameter, 
S(k, t) = (ipzipz), averaged over all possible orientations 

of k and 50 independent initial conditions. Such study is 
standard If \k - fc | < fc and k a L(t) > 1, S(k,t) 

satisfies the scaling form S(k,t) = L(t) f[(k 2 — 
k 2 )L(t)Xo], where d the spatial dimension. Analysis of 
the moments of S(k,t) shows that L(t) ~ t 1 ^ with 
1/z ~ 0.31. We also find that S(k ,t) ~ t - 32 . 

Finally, we have verified that the value of the exponent 
z calculated from cither the grain boundary perimeter or 
the structure factor is not modified by the introduction 
of random fluctuations into Eq. ([!]) . 



The results presented are qualitatively modified fur- 
ther from onset. For example, at e = 0.25 we find that 
different linear scales of the structure are no longer pro- 
portional to each other, and we obtain effective exponents 
that are in agreement with the values of z reported in ear- 
lier studies at e ~ 0.25 ||J|. We find 1/z ~ 0.21 from 
an analysis of the moments of S(k,t), 1/z ~ 0.26 from 
the grain boundary perimeter, while the moments of the 
distribution of curvatures yield 1/z ~ 0.32. 

We next discuss a possible growth mechanism lead- 
ing to an exponent of z = 3, as well as our interpreta- 
tion for the slower growth that is found when e is not 
sufficiently small, and its dependence on random fluc- 
tuations. Coarsening exponents can be often inferred 
from the law of motion of the class of defects that con- 
trol coarsening |19|. A typical transient configuration 



(Fig. [I]) contains a large amount of grain boundaries, as 
well as other defects such as +1/2 disclinations. Grain 
boundaries move over large distances, whereas disclina- 
tions remain largely immobile. Disclinations produce a 
background of curved rolls that cannot freely relax due 
to topological constraints (this is in contrast with other 
coarsening mechanisms discussed in ^,0].) However, at 
any given point roll relaxation as well as disclination an- 
nihilation do occur after the passage of a grain boundary. 
In ref. p5| we studied the motion of a grain boundary sep- 
arating two semi-infinite domains of mutually perpendic- 
ular rolls, straight on one side, curved on the other. At 
lowest order in e, those parallel to the boundary are dis- 
torted whereas those that are perpendicular to it remain 
straight. The energy of the configuration decreases by 
a net displacement of the grain boundary, the effect of 
which is to replace curved rolls by straight ones of lower 
energy. Therefore the size of the domain with curved 
rolls decreases. We showed that if the curvature of the 
rolls ahead of the boundary is k, the boundary advances 
at a speed, 



v gb 



(3) 



It was shown in @ that Eq. (|) is in quantitative agree- 
ment with a direct numerical solution of Eq. (|l|) with an 
initial condition that involves a 90° grain boundary. As 
seen in Fig. [I], disclinations produce roughly axisym- 
metric patterns, with a characteristic roll curvature that 
is inversely proportional to the distance among disclina- 
tions. If this distance is proportional to the characteris- 
tic linear scale R(t), then dimensional analysis of Eq. (||) 
suggests R(t) ~ t 1 / 3 . 

These considerations are modified further from onset. 
The law of grain boundary motion, Eq. (^), is valid only 
to first order in e. At moderate values of e, the sepa- 
ration of length scales assumed in the derivation of Eq. 
(^) breaks down (the grain boundary thickness is of order 
Xo/y/e), leading to non-adiabatic effects. Within the am- 
plitude equation formalism, and following the approach 



3 



of ref . [gCj , we have obtained the leading order nonadia- 
batic corrections, and find that the position of the grain 
boundary and the phase of the rolls located ahead of it 
no longer decouple. Equation (p|) can be generalized to, 



v gb 
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D(e) 
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where Xb g is the average location of the grain boundary, 
(/> is a constant phase and, 
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p[e) = maxe 



dx AQ(x)d x A (x) cos(2fc a; + 9) 
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+ - / [2A QJ Bga x ^ + A 2 B d x B ] cos(2fc a; 
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The functions Aq (x) and Bq (x) are the amplitudes of the 
two sets of rolls separating a planar boundary the 
coefficient -D(e) represents a friction term, and p(e) the 
amplitude of a periodic pinning potential. The contribu- 
tion from non adiabatic effects is typically of the order 
of, 
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where \a\ is a constant of order unity. Hence, p behaves 
non analytically as e — > 0, but increases extremely fast 
with e at low values of e. From Eq. (||) we see that 
at any finite e > 0, there exists a critical curvature n g 
below which v g b = 0. Remarkably, pinning becomes no- 
ticeable even at e = 0.1 p5| , and grain boundaries were 
seen to advance only by half-integer multiples of the roll 
wavelength. Therefore we expect that grain boundaries 
in a coarsening configuration will become pinned over 
time. We believe that this pinning is the reason behind 
the lower effective exponents found in previous studies 
at e = 0.25 [^|joK^], as well as for the related result that 
random fluctuations added to Eq. ([[]) consistently lead 
to larger coarsening rates. 

In summary, we have presented results for several in- 
dependent measures of the linear scale of stripe patterns 
ordering very near onset, and obtained a coarsening ex- 
ponent that is very close to z — 3. This value can be 
explained through dimensional analysis of the velocity of 
a single grain boundary advancing into a background of 
curved stripes in the limit e — > 0. This mechanism also 
predicts increasing corrections to scaling further from on- 
set due to defect pinning. 
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